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Ch ■ Abstract 

O 

J^ , We consider wave propagation inside an anisotropic superconducting film 

m . 

f^ , sandwiched between two semi-infinite non-conducting bounding dieletric me- 

(N 

^^ \ dia such that along the c-axis, perpendicular to the surfaces, there is a plasma 

frequency Wp below the superconducting gap. Propagation is assumed to be 

parallel to the surfaces in the dielectric medium, where amplitudes decay ex- 



'^ \ ponentially. Below ujp, the amplitude also evanesces inside the film, and we re- 



trieve the experimentally measured lower dispersion relation branch, lo oc v^) 
and the recently proposed higher frequency branch, to oc l/\/^- Above Wp, 



o 
o 

> 

^. 

H , propagation is of the guided wave type, i.e., a dispersive plane wave confined 

inside the film that reflects into the dielectric interfaces, and the modes are 
approximately described by a; « Up^Jl + (/3//3o)^, where /3o is discussed here. 



I. INTRODUCTION 

The experimental measurement of a plasma edge in the infrared reflectivity of High- 
Tc bulk superconductorsETu La2-xSrxCu04, YBa2Cu30s-xO, and Bi2Sr2CaCu20^U has 
brought a renewed interest in collective oscillations with the Cooper pair density. Previ- 
ously, such studies were hampered by the well-known arguments that the Coulomb interac- 
tion shifts the frequency of such oscilations to above the gap frequency. Recent theoretical 
studiesli3^tJ support the view of a plasma oscillations along the c-axis, the direction orthog- 
onal to the Cu02 planes, caused by a large superconducting gap and a high anisotropy in 
these materials. 

In superconducting films, plasma oscillations have distinct properties. They exist in films 
regardless of their anisotropy or layered structure, and regardless of the type of pairing or 
even of the critical temperature value. This mode was predictedll^'tZlll^ some time ago and 
was bothly observed in thin granular aluminium films, in the hundreds of MHz rangecJ, 
and in thin YBa2Cu^0s-x filmstj, in the higher frequency range of hundreds of GHz. 

The study of propagating modes in films is not a new research field, they have been 
measured since long ago, in several materials ranging from metalsc3c2l to semiconductorsc3, 
and also discussed in the theoretical literaturcMMEa. The novelty for the superconducting 
film is the very low frequency where these modes are observed and the strong temperature 
dependence, explained by a collective oscillation with the Cooper pair density. RecentlycZl it 
has been proposed that, similarly to metals and semiconductors, two branches of propagat- 
ing modes should be observed in highly anisotropic superconducting films, e.g., the High-Tc 
compounds. The experimental observation of these two modes can provide a method to 
independently measure the London penetration lengths, parallel and perpendicular to the 
film surface. 

In this paper we include into the study of mode propagation in films the existence of a 
c-axis plasma frequency, tUp, inside the superconducting state. We only consider the c-axis 
perpendicular to the surfaces since this configuration is the most easily grown nowadays. 



The propagation of plasma modes above and below uop, their dispersion relations are our 
goals in the present paper. We restrict our study to waves that propagate along the surfaces 
of the superconducting film and evanesce in the dielectric medium perpendicularly to the 
surfaces. Two very distinct behaviors are predicted, namely, field amplitudes that also 
evanesce inside the film away from the interfaces {uo < Up), and confined propagation, i.e., 
waves propagate inside the superconductor along an oblique direction and reflect at the 
dielectric-superconductor interfaces {lu > UpfB. In many aspects this last case resembles 
propagation in optic fibers. In conclusion this paper considers a plasma frequency Up inside 
the superconducting state and studies its effects into wave propagation in films. 

The study of plasma modes in superconducting films with a plasma frequency Up has 
been previously considered by Artemenko and Kobel'kovEa in the context of kinetic equations 
for Green functions generalized to the case of layered superconductors with weak interlayer 
coupling. The present paper provides a more complete study of such modes than theirs 
because it shows that there are symmetric and antisymmetric modes for uj < Up and dis- 
cusses their several distinct regimes. The antisymmetric mode is intimately connected to 
the transverse current component, thus being highly sensitive to the transverse London pen- 
etration length, in the same fashion that the symmetric mode depends on the longitudinal 
London penetration length. We claim that these properties can be used to gain a better 
understanding of the transverse and longitudinal supercurrent components in anisotropic 
and layered films. Besides, the present study is done in a framework distinct from those 
authors, we take the simplest possible theory, the London-Maxwell theory for anisotropic 
materials, supposedly valid when all lengths are larger than the inter-layer separation a. 
In the context of this theory we determine exactly the mode frequency, u;, as a function 
of its wavelength along the surfaces, (3, and of the remaining parameters: the two London 
penetration lengths, transverse (A_l) and longitudinal (Ay) to the surfaces; the dielectric 
constant of the non-conducting medium exterior to the film, e; the film thickness d; and 
the isotropy and frequency independent Eg, the simplest phenomenological choice for the 
superconductor's static dieletric constant. 



This paper is organized as follows. In the next section, (||), we derive the dispersion 
relation equations, above and below the plasma frequency, using the London- Maxwell the- 
ory, and solve them exactly. Our study is restricted to identical top and bottom dielectric 
media. Similar conclusions should also apply to the most general asymmetric case. Section 
pn] deals with the u < Up case, where we study the two possible branches, the so-called 
symmetric (lower) and antisymmetric (upper) branches, and their three possible regimes: 
optical, coupled, and asymptotic. For each branch we derive, from our exact solution, useful 
approximated expressions for each of the above regimes. In this fashion many of the previous 



film studies are shown to be suitable approximations of our exact solutions. In section ^ 
we study wave propagation above Up and, obtain, from our exact dispersion relation, the 
several dispersion relation branches found by Artemenko and Kobel'kovE^. We show that 
symmetry still plays an important role and the two branches below Up split into several 
ones, the symmetry of the state being determined by the number of half-wavelengths that 
fit perpendicularly to the film. In section |V| we apply the present theory to the High-Tc 
superconductors, suitably choosing the parameter values. Finally in section |VI| we summa- 
rize our major results. The proof of some of results from our exact solution, such as the 
approximated expressions, are the subject of appendix |V11. 



II. THE LONDON-MAXWELL THEORY APPLIED TO ANISOTROPIC 

SUPERCONDUCTING FILMS 

In this section we apply the London-Maxwell theory to describe wave propagation in a 
superconducting film sandwiched between two identical non-conducting dielectric mediaa. 
An external electromagnetic wave of frequency u, and wavenumber k = uj/c, is absorbed by 
the superconducting film, producing a mode whose dispersion relation is uj{P). We choose 
a coordinate system where the two plane parallel surfaces separating the superconducting 
film to the dielectric medium are at x = d/2 and x = —d/2 . Propagation takes place along 
the z axis such that all fields can be expressed as Fi{x) exp [—i{Pz — ut)]. At this point 



we introduce the time dependence exp(i uj t) to all fields into the basic equations governing 
the system. Current transport inside the film is described by the first London equation, 

i a; /ioA| J|| = E||, i u noX\ J± = E±, (1) 

where Ey and E^ are the field components parallel and perpendicular to the film surfaces, 
respectively. The electromagnetic coupling, given by Maxwell's equations, 

V ■ D = (2) 

V- H = (3) 

V X E = -iujfioH (4) 

V X H = iuD where D = e^E — iJ/u. (5) 

shows that the superconductor dielectric constant is tensorial, D = eQE ■ E. 
fe± 0\ 

ei\ 
VO eii/ 

Phenomenological theories, such as the present, describe the superconductor only at 
energies much lower than the pair breaking threshold. The frequency u is much smaller 
than the frequency defined by the superconducting gap. Consequently the frequency is 
much smaller than the plasma frequency along the Cu02 planes, k\\\ <^ 1, thus rendering a 
negative, and large in modulus, dielectric tensor component parallel to the surfaces. 

In the present model the c-axis plasma frequency is reached when the dielectric tensor 
component perpendicular to the surfaces becomes null. 

e±{uj = ujp) = LUp = ——— (8) 

Hence e^ changes sign as u crosses Up. 
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eL = es- JT^^^ £\\ = £s - JZV^^ (6) 



Solving Maxwell's equations for this particular geometry gives two independent sets of 
field components, {Hx,Ey,Hz), the transverse electric (TE), and [E^^Hy^Ez)., the trans- 
verse magnetic (TM) propagating modes. For the TM mode the interfaces superconductor- 
dielectric acquire superficial charge densities. Curiously anisotropy also implies in a volu- 
metric charge density inside the film, p = — iV ■ J/u = eo(l — e±/e\\)dE±/dx, although it 
is not responsible for these low frequency propagating modes. In fact such modes were first 
observed in nearly isotropic thin granular aluminum filmsta. In conclusion the TM mode 
supports low frequency travelling waves because of the oscillating superficial charge densities 
that couple the non-conducting dielectric medium to the superconducting film, their ratio 
given by, 

a{x = d/2) _ dEjdx\^=d/2 /qv 

aix = -d/2) ~ dEjdxU=.,/2 ^ ^ 

The TM field equations are given below, for the dielectric medium, {x > d/2 and x < —d/2), 

T^ OX T^ OX OX^ 

and for the superconducting strip {—d/2 < x < d/2), 

^ 3e\\ dEz ,, . uJEw dEz d'^E^ ,„ „ 9 £11 ^9 ,9 ,,,, 

T^e± ox r^ ox ox^ e± 

In this paper we seek propagating modes that evanesce in the dielectric medium, namely 

display exponential decay in the dielectric. This condition is f^ = /?^[1 — (uj/Pv)'^] > 0, and 

can be interpreted as demanding a phase velocity, uj/jS, smaller than the speed of light in 

the dielectric, v = c/ \[e. Hence above the film [x > d/2) one gets that, 

Ez = Kexp{-fx) (12) 

From the other side exponential decay inside the film is not sure to occur because the sign 
of r^ is not uniquely defined. This sign determines different physical regimes, and for this 
reason, we introduce the below and above plasma condition and study their respective r^ 
regions. Hereafter, in order to simplify further discussions, we will choose the nonconducting 
dielectric medium to have the largest static constant, e > Eg- 
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UJ < 00- 



p 



In this case both e\\ and e± have the same sign. Introducing the Eq.(^ ap- 



proximation, one gets that r^ ^ k||/£^±|/3^ + iZ-^y > 0- The planes of constant phase are 
z = const, and field amplitudes also evanesce inside the superconductor. 



^. 



Eg exp (— r x) + Fo exp (r x), r^ > 



(13) 



UJ > Ur, 



In this case e\\ < and e^ > 0. Even within the Eq.(|^ approximation, 



r^ ^ — k||/^±|/3^ + 1/Am has no definite sign. The mode may be evanescent or propagative 
inside the superconductor. However the evanescent modes disappear for e > Eg and then 
one gets that r^ < 0. In this case the planes of constant phase are /3z ± r'x = const. 
{t'"^ = — r^ > 0), and we fall into the case of confined propagation, with an oblique incidence 
well defined at the interfaces. 

£11 



Ez = Eoexp {-iT x) + Fo exp {ir' x), r ^ = -^(3^ + Pe\\ > (14) 

In summary for the present purposes the sign of r^ is uniquely determined below and above 
Up by the sign of the ratio e\\/e±. 

The dispersion relations follow from the continuity of the ratio Hy/Ez at a single interface, 
say X = d/2, once assumed that the longitudinal field Ez is either an even or an odd function 
with respect to the x = plane. This is possible because the two dielectric media, above 
and below the film, are equal. 

Table I : The four possible dispersion relations and Ez. 



Ez 


Dispersion Relation 


Eoz cosh (r x) 
Eoz sinh (r x) 
Eoz cos (r' x) 
Eoz sin (r' x) 


ff- = -tanh(r f) 

re _ 1 


r £|| tanh(r |) 

f^ = -tan(r'f) 

t' e _ 1 


f e\\ tan{r' |) 



(15) 

At this point we introduce some dimensionless variables useful in the study of the above 
dispersion relations. 



Table II : Dimensionless Variables. 



Retardation 



7 



Ujjfi 



Anisotropy 



R = {^? 



Dielectric 



£3 



Thickness 



^-(^)' 



Wave Number 



x = (/?A|,r 



(16) 



We cast some of our previous results in the dimensionless variables formalism. The condition 
of evanescence along the direction orthogonal to the surfaces is < 7 < 1, since according 
to Eq.(E 



, one has that E^ = Egexp (— y ^ a/1 — ^y'^x/d). The dominance of the incident 
wavelength over London's penetration along the surface, which gave Eq.(^, becomes 7^ -C 
r/X. The dielectric component ratio, within the Eq.(|^) approximation, is b±/6\\ ~ (1/-R) (1 — 
i?7^ X/r). Thus the u < Up and u > Up regimes correspond to R'y'^ X/r < 1 and i?7^ X/r > 
1, respectively. The c-axis (High-Tc compounds) is perpendicular to the surfaces and so 
R>1. 

In the following, we obtain the exact solution for each of the Eq. (|I5|) dispersion relations 
considering the dimensionless variables dr {uj < ujp) and dr' (uo > uOp) as a parameter s, 
whose range is still to be determined below and above ujp. Therefore we seek the dimension- 
less curve, 7(X), in parametric form, [X(s),7(s)], and immediately obtain the dispersion 
relation. 



{s)Ms)] = \4x{s)/\,{v/\i)^{s)^)]. 



uj < ujp We introduce the dimensionless parameter. 



t = dr 



XR 



+ 1 



Using the Eq.(^) dimensionless variables, the Eq.(|T5|) dispersion relations become, 



(17) 



Table III : The dispersion relations for uj < ujp 



E._ 



Eoz cosh (tx/d) 
Eoz sinh (tx/d) 



Dispersion Relation 



VAX 



v^ 



tanh {t/2) 
t 

1 



V^l^ 



t tanh (t/2) 



(18) 

To find their solution, first obtain X(t, 7) from Eq.(|l^, and then introduce it into Eq.(|l^). 
We obtain a second degree equation for 7(t)^ with two roots: one negative the other positive. 
The negative root does not correspond to propagating modes and only the positive is left. 



7/(t) 



^ (^At''-l)-r + ^[{At^ _ 1) + r]2 + 4 ^ (At2 - 1) gj{t) 



2{At^ 



Dfl + 



Ar 



9iit)] 



X(t) 



At^-1 



(19) 
(20) 



/ = {S,A} labels the E^ symmetry: 5f^(t) = t^ tanh^ (t/2) and gsit) = tV^anh^ (t/2) are 
associated to the antisymmetric and the symmetric modes, respectively. According to the 
above equations the parameter t range is [l/vA, cxd]. 



u > Up Through Eq. ([I^ we define a new dimensionless parameter. 

1 



t' = dT' 



^\ 



XR 



1 



(21) 



r 

Introducing the Eq.(0) dimensionless variables into the Eq . (pTsD dispersion relations gives 
that, 

Table IV : The dispersion relations uj > ujp 



E,, 



Eoz cos {t'x/d) 
Eoz sin (t'x/d) 



Dispersion Relation 



VAX^ 



Vi^ 



tan(i72) 
t' 



VAX 



1 



v^ 



t' tan (472) 



(22) 



To solve them, express Eq.(0) as X(t',7) and introduce it back into Eq.(p2|) obtaining a 
second degree equation for ■y{t'Y, whose solution is. 



li{t 



Xiit') 



{At'^ + 1) + r ± V^[(v4t'2 + 1) - r]2 - 4 ^ [At'^ + 1) hi{t') 



2(At'2 + l)[l + ^hiit')] 



At'^ + l 



R-ii{t'f (At'2 + l)-r 



(23) 
(24) 



Like in the previous case, I = {S,A} gives the Ez parity: hA{t') = f^ tan^ (t'/2) and 
hs{t') = t'^/tan^ (^V2) ior the antisymmetric and symmetric modes, respectively. 

The discussion of the suitable parameter range is more involving in the present case. 
According to Eq.p^) the function t'tan(t'/2) must remain negative for the antisymmetric 
problem whereas tan(t'/2)/t' must be positive for the symmetric case. This restricts the 
parameter range t' to the intervals, [2 A^ vr, (2 A^+ 1) vr] for symmetric and [(2 A^+ 1) vr, 2 (A^ + 
1) tt] for antisymmetric, where A^ is an integer larger or equal to zero. 

Our exact parametrized solutions uj{(3) are ready for applications, just requiring the 
numerical parameters. Before doing so in section 0, we find convenient to describe their 
physical properties in the next two sections. 

III. DISPERSION RELATION BELOW THE PLASMA FREQUENCY 

In this section we study the properties of propagating modes below ujp. We find con- 
venient to summarize their major physical properties and introduce some approximated 
expressions, each describing a different regime of the exact uj < ujp dispersion curve. We 
leave to the appendix [VII| the proof that the previous section exact results do justify our 
picture and yield the approximated expressions. Below the plasma frequency the field ampli- 
tudes inside the film evanesce from the surfaces, and this exponential decay is characterized 
by r, according to Eq.(|l^). The choice of film thickness, d, is very important in order to 
assure sufficient coupling between the two surfaces. For extremely thick films the surfaces 
decouple and the symmetric and the antisymmetric modes are the same, only independent 
surface plasma modes exist in this situation. Indeed, in case of strong coupling between 
the two surfaces three distinct regimes are possible for both symmetric and antisymmet- 
ric modes. In sequence of increasing /3 we call them optical, coupled and asymptotic (see 

Fig.(l)). 

Close to the origin /? ~ 0, the modes are optical, that is, they are essentially plane waves 
travelling in a dielectric medium. There is almost no exponential decay in the surrounding 
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dielectric semi-spaces up to appreciable distances, f ^ 0, according to Eg. (p!2D . Propaga- 
tion occurs with the speed of light in the dielectric, (7 ~ 1). The kinetic energy of the 
superconducting carriers is negligible compared to the magnetic energy of the mode. The 
superconducting film contributes very weakly to this regime. In the optical regime one has, 
approximately, the linear behavior u ^ v (3. 

Slowly the condensate's kinetic energy increases with (3 until finally both energies become 
comparable and a first cross-over takes place. This is the onset of the coupled regime, a 
slow mode (7 < 1), where the condensate's kinetic energy dominates over the magnetic 
energy, the film and the dielectric are strongly coupled. This is the true plasma mode 
regime. Such evanescent propagating modes are knownEJllZlEJ, here we discuss the effects 
of a c-axis plasma frequency Up below the gapBa. In the coupled regime inside the film, the 
fields evanesce from the surfaces very smoothly resulting into two types of coupling between 
the superficial plasmons, which are the antisymmetric and the symmetric modes. For the 
symmetric mode its dispersion relation has been theoretically studied by many authorsEjllH 
and experimentally observedli^ilHl in the past: 



V IdB ,^^. 

Recentlynll it has been proposed that the symmetric mode is just the lowest frequency branch. 
In fact, there is an upper (antisymmetric) branch mode, that can also be experimentally 
observed for highly anisotropic High-Tc superconducting films: 



Because the antisymmetric branch is the highest in frequency, it is more sensitive to the 
effects of the plasma frequency Up. For di3/2 comparable (or smaller) to Ss/e, this branch 
becomes difficult to observe since uja ~ ujp. For dj3/2 dominant over Ss/e, the symmetric 
and antisymmetric relations are proportional to a/5 and l/y/P, respectively. 

The cross-over frequency between the optical and the symmetric regime is obtained at 
the intersection of the two branches: 
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vd 
2A1 



^cross,s ~ o\2 V / 



The cross-over between optical and antisymmetric modes is implicitly given by the unique 
real solution: 

(^) ^lross,a + (j) ^loss^a ' i^f = (28) 

For the antisymmetric mode, the optical and the coupled frequencies are, respectively, an 
increasing and a decreasing function of /3. Thus the cross-over uJcross,a also gives an estimate 
of the maximum frequency, associated to the peak seen in the antisymmetric curve (Fig.(^). 
By increasing (3, evanescence inside the superconductor becomes stronger (r is large), up 
to the point where the surfaces are nearly decoupled. We have just reached the asymptotic 
regime. Thus for sufficiently high /3, symmetric and antisymmetric modes converge to the 
same asymptotic frequency. 



f^as — ^p 



l + Jl + A^'"'^' 



\ l + 2(|^)^ + Jl + 4(--^ 



) (29) 



Notice that it is always true that uJas < ujp. It is easy to check in the following extreme cases 
of a very high dieletric constant {e/Es ^ ^±/^\\ )> and of the opposite case, namely, when 
anisotropy plays a more important role than the dielectric constant, {e/Es <^ '^±/'^|| )• In 
the former case we obtaintH uJas = v/J\±\\\, and in the latter cUas = ^p- In summary we 
have just reviewed the major features of the three possible regimes for both the symmetric 
and the antisymmetric dispersion relations. 

We find that the dielectric constant ratio r must be sufficiently high in order to prevent 
that the (3 range be extremely limited for the coupled regirae. In the next section (0) we study 
such crossovers in case of the High-Tc parameter values. The criterion for the disapperance of 
the coupled antisymmetric and symmetric regimes is just given by an asymptotic frequency 
Uas smaller than ujcross,s and uJcross,a, respectively. As an example consider the case of an 
extremely large anisotropy (Aj_ — » oo), where we get that Uas ~ cUp ^ 0, eventually smaller 
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than the crossover frequencies of Eq.(p7D and Eq.(]28|). In this hmit the dispersion relation 
goes almost directly from the optical to the asymptotic regime. 

IV. DISPERSION RELATION ABOVE THE PLASMA FREQUENCY 

In this section propagating modes inside the superconducting film with the frequency 
larger than ujp are considered. We summarize their major physical properties, and leave 



to the appendix [VI I| their derivation from the exact parametrized dispersion relation. We 
are again interested in modes that evanesce in the dieletric medium. As discussed before, 
r' is always real with the experimental parameter choice of e > Eg- The propagation of 
light inside the film can be understood by a superposition of two dispersive plane waves, 
exp {— 'i[±r'(/9)x + [3z — co'(/?)t]}, regarded as incident and reflected waves, thus with a well 
defined angle of incidence at the interfaces superconductor-dielectric. Thus while outside 
the film light is still evanescent, in its interior the propagation can be pictured through ray 
optics (see Fig.(|I])), similarly to an optical fiber. 

The exact parametrized plasma relation, given by Eq.(^) and Eq.(^), leads to the 
simple approximated dispersion relation, in the vicinity of uj/ (3 ~ w. 



UJm — ^p 



{/3dy A^ 2 

where M is a positive integer, odd for the symmetric modes (M = 2N + 1) and even for the 
antisymmetric modes (M = 2(A^ + 1)). This approximated expression was first obtained by 
Artemenko and Kobel'kovE^. The M*'* curve satisfies the evanescence condition, 7 < 1, for 

l3init{M)= , (31) 



e (7rM)2VA||/ 

At the particular value (3 = (3init{M), the wave vector along the x-direction is given by 
t' d = Mtt where M is exactly the number of half- wavelength that fit perpendicularly to 
the film. For j3 > l3init{M), the relation between r' and M is not so simple because the wave 
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is also in the dielectric medium. In this case M just determines the number of E^ extrema 
along the x-direction (see Fig.(|I])). For instance, when just one half- wavelength fits into 
the film we are facing the M = 1 symmetric mode, which has just a single Ez extremum 
(maximum). We notice that the frequency of the plasma modes decreases for increasing M. 
For (3 < (3init{M) we have that 7 > 1 and the modes are not evanescent, they are plane 
waves travelling in the dielectric medium. They have an oblique incidence at the surfaces 
and to see this just take the plane wave exp {—i[fx + j3z — uot]} with the speed of light in 
the dielectric, v = u/y/f^ + /3^. One obtains that f = /3-\/7^ ~ 1; showing that we are in a 
7 regime not studied here. 

The present theory eventually breaks down because large r' means an infinitesimally 
small wavelength inside the film along the x axis. For this reason a physical cut-off must be 
put into this theory, the inter-plane separation a. The upper limit r/^„^ is determined when 
a x-direction half-wavelength fits into the interplane distance a. 

'^Lax = 7r/a = M^axT^/d (32) 

Introducing M^ax in the approximated dispersion relation Eq. (PP]) , gives the lowest fre- 



quency plasma mode: uJmin = i^p y 1 + (/^Aj/Tr)^, where \j = aAj_/A|| is the Josephson 
penetration depth. 

Beyond t^^x niore elaborate models should provide a description of the condensate in 
such scale. This can render the theory's usage quite limited, in case the ratio d/a is small, 
because for very few half- wavelengths inside the film we reach the cut-off limit. However 
the applicability of this model above Up is not restricted to a small ratio d/a because the 
propagating mode has constant amplitude inside the film, and so, no matter how far apart 
the surfaces are, there will always be waves reflected at the surfaces. For this reason the 
major conclusions of this section, namely, the existence of slow modes above the plasma 
frequency, corresponding to dispersive plane waves inside the film, must remain valid even 
near the cut-off limit. 

In the next section we take the standard parameter values for the high-Tc material and 
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discuss the properties of such propagating modes below and above the plasma frequency. 

V. APPLICATIONS 

For our applications we choose the following set of parameters for the High-Tc ceramic 
superconductors. The static dielectric coustauti is .. « 30. The zero-temperature London 

penetration length along the Cu02 planes is Xco = 0.15/i?7TEj. While the anisotropy, A_l/A||, 
is 5 for YBa2Cu30s-J^, it has been changing in the past for Bi2Sr2CaCu20x, ranging 
from 10^, mainly from torque measurements studie£ll, to much higher valuesE3. We show 
here that choosing the anisotropy between these two values can have important effects 
on the properties of the dispersion relations for uj < ujp. There is also the compound 
Tl2Ba2CaCuO^ compoundU with \^/\\\ ^ 90. 

Fig.(|l]) provides a pictorial intuitive view of wave propagation in the film below and above 
the plasma frequency. Below ujp the instantaneous electric field, as well as the component 
Ez, are shown for both symmetric and antisymmetric fields. The surface charges are the 
sole sources of propagating electric fields. The superficial charge arrangement has strong 
consequences for the electric field distribuiton inside the film, leading to symmetric and 
antisymmetric modes, found at distinct frequency ranges, the latter being an upper mode. 
To understand the effects of anisotropy, or a layered structure, into these modes, where the 
Cu02 planes are parallel to the film surface, consider the transverse and longitudinal field 
and supercurrent components. E^ and Jx are very intense for the antisymmetric mode and 
nearly zero for the symmetric one, and E^ and Jz-, are the dominant components for the 
symmetric mode but not for the antisymmetric one. Thus the relevant penetration depths 
for the symmetric and the antisymmetric modes must be Ay and A_l, respectively. The 
anisotropy A^ > X\\ hardens the system along the c-axis thus making the antisymmetric 
mode lower in frequency. Above Up a dispersive plane waves travels inside the film that 
undergoes total reflection at the interfaces. This Figure also depicts the number of extrema 
for the electric component E^, which determines the symmetry of the mode. 
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In Fig.(|^), the symmetric and antisymmetric dispersion relations are shown for a very 
thin film d = lOnm and the three anisotropies, X±/X\\ = 5, 10^ and 10^. In order to 
slow down as much as possible light in the dielectric, and so, lower the coupled regime 
frequency range, we choose SrTiOs as the exterior non-conducting dielectric medium. At 
low temperatures its dielectric constant is known to be high up to the GHz frequencytS, 
e ~ 20000. All curves in this Figure were directly obtained from the parametrized solution of 
Eq.(p!9l) and Eq.(|^). The symmetric state crossover frequency, u!cross,s, between the optical 
and the square root regimes is found to be 75 GHz for the three anisotropies displayed 
here. The antisymmetric state crossover frequency, u>cross,a, between the optical and the 
coupled regime, depends on the anisotropy, being 2.3 THz, 290 GHz and 48 GHz for 
'^±/'^|| = 5, 10^ and 10^, respectively. Because of this anisotropy dependence they are not 
indicated in this Figure. Notice that the symmetric mode (dashed lines) is always found at a 
frequency range lower than its corresponding antisymmetric mode (continuous lines). In the 
low wave vector limit, all the dispersions collapse into the same linear (optical) regime. In 
the opposite limit, (3 very large, both symmetric and antisymmetric curves converge to the 
same asymptotic frequency. This asymptotic frequency is just the surface plasma frequency, 
strongly affected by anisotropy: 1 THz for YBCO, 200 GHz for BSSCO and 50 GHz 
for A_l/A|| = 10'^. This saturation frequency is independent on film thickness, as expected, 
according to Eq.(|29D. Notice that the wave vector signaling the onset of the asymptotic 
regime decreases for increasing anisotropy, ranging from lOfim'^ for YBCO to lOOmm"^ for 
the maximum anisotropy considered here. Between the optical and the asymptotic is the 
coupled regime, the true plasma mode, found to extend over a large range of frequency and 
wavenumber for an YBCO thin film, according to this Figure. In this regime the symmetric 
plasma mode follows a square root and the antisymmetric mode clearly shows the inverse 
square root dependence. For BSCCO, the coupled regime is shorter than in YBCO, ranging 
between 0.1 to 5fim~^. For the maximum anisotropy of 10^, the coupled regime disappears 
thus only surviving the optical and the asymptotic regimes. We conclude that an extremely 
large anisotropy inhibits coupling between the two surfaces and enhances the surface plasma 
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modes. 

Fig.(H), directly obtained from Eg. (pQ]) , shows the asymptotic frequency versus e, for 
several anisotropies. Notice that the asymptotic frequency for small e is just Up. For 
YBCO, the asymptotic frequency drops over an order of magnitude when e changes by 
three decades. The choice of a non-conducting medium of small dielectric constant can 
render impossible the observation of these plasma modes because the asymptotic frequency 
becomes larger than the superconducting gap. In BSCCO, and other extremely anisotropy 
compounds, the asymptotic regime is always much smaller than the gap. As e increases, we 
have determined the crossover value, esX±/X\\, where the assymptotic frequency acquires a 
e dependence, thus being strongly affected by anisotropy. Hence in order that the inverse 
square- root dependence of the antisymmetric coupled regime be observable, the condition 
e > esX±/\\\ must be satisfied. 

Fig.(^) shows the dispersion relations above Up for a 100n?Ti-thick film. The dielectric 
constant e is taken equal to Es- We have chosen a BSCCO compound with an anisotropy 
'^±/'^|| = 10^- Therefore the c-axis plasma frequency is much smaller than the superconduct- 
ing gapEj. The curves were directly obtained from Eq.(PBD and Eq. (P^) . Above cUp, plasma 
modes branch into a large number of modes when film thickness increases. For the choice of 
a 100n?7i-thick film, the present London-Maxwell theory remains valid for the first 66 modes. 
We plot in this figure the M = 1, 2, 3, 4, 7, 12 and 66 dispersion relations. For increasing M 
the plasma mode becomes slower, indicating that the number of reflections that the confined 
plane wave undergoes at the interfaces has also increased. Notice that the M = 66 plasma 
mode is very near Up. The M=l branch is very near to the optical branch showing that it 
essentially has the speed of light in the dielectric medium. 

VI. CONCLUSION 

In this paper we have studied propagating plasma modes in a superconducting film 
surrounded by two identical dielectric media. The superconductor is anisotropic having its 
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uniaxial direction (c-axis) perpendicular to the interfaces with the dielectric medium. We 
consider the existence of a plasma frequency along the uniaxial direction below the gap and 
study its effects into the propagating modes using the London-Maxwell theory, which gives 
a good account of the physical situation for scales larger than the interplane separation. 
In fact we do use the interplane separation as a cut-off for the present theory. We only 
consider low incident frequency compare to the superconducting gap frequency. Moreover 
the wavelength associated to the incident wave is much larger than London's penetration 
length along the surfaces. Therefore the dielectric constant of the superconducting film along 
the plane is always negative and with a large modulus. Along the c-axis, the situation is quite 
distinct. The dielectric constant of the superconducting film along this direction vanishes 
at the plasma frequency, being positive above. In this paper we were only concerned about 
modes that propagate parallel to the surfaces and evanesce in the dielectric medium. Under 
this condition, and the requirement of the non-conducting medium dielectric constant larger 
than the superconductor's static constant, we find that modes are either totally above the 
plasma frequency, or below, thus never crossing the plasma frequency line. Thus we study 
these two regions of frequency separately. Within our London-Maxwell theory framework 
we find the exact expressions for the dispersion relations of the plasma modes in the two 
cases. The exact expressions are found in parametric form, the parameter being the one that 
characterizes the exponential behavior inside the superconductor. Approximated expressions 
for some especial regimes of the dispersion relations are obtained from our exact parametric 
solution. In this way we retrieve well-known results in the literaturecJilllOEa to the context 
of an anisotropic film with a plasma frequency inside, and also, derive new ones. We find 
that propagating modes below and above the plasma frequency have quite distinct physical 
properties. 

Below the plasma frequency, the amplitude also evanesce inside the superconductor, the 
film thickness is important in order to assure a sufficiently strong coupling between the sur- 
faces. There are two branches of the dispersion relation corresponding to the two possible 
arrangements of the superficial charge densities. The symmetric branch is the lowest in 



frequency and has the two superficial charge density symmetrically disposed. The highest 
branch has opposite charge facing each other at the interfaces thus being antisymmetric. We 
have studied here in details the three possible regimes that can exist for these two branches, 
their cross-over, and the conditions for the existence of the coupled regime, the most interest- 
ing of the three regimes. In the so-called coupled regime both symmetric and antisymmetric 
dispersion relation branch provides independent information on the two London penetration 
lengths. The antisymmetric mode is intimately connected to the transverse current compo- 
nent, thus being highly sensitive to the transverse London penetration length, in the same 
fashion that the symmetric mode depends on the longitudinal London penetration length. 
We believe that this remarkable property can be used to gain further understanding of the 
transverse current component in anisotropic and layered films. 

Above the plasma frequency the amplitudes do not evanesce inside the film, propagation 
is geometrically understood, and an oblique ray can represent the plane wave inside the 
superconductor which is totally refiected at the interfaces (see Fig.(|I|)). However outside the 
film, in the dielectric medium, the wave remains evanescent, and we conclude that above 
the plasma frequency the superconducting film displays confined propagation, similar to an 
optic fiber. Notice that the film thickness is not a crucial parameter because no matter how 
apart the surfaces are, there will always be plane waves travelling in its interior. Contrary 
to below the plasma frequency, above there are many modes corresponding essentially to 
the number of half-wavelengths that fit inside the film along the c-axis. 

In summary we have shown in this paper that superconducting films surrounded by a 
dielectric medium displays very interesting plasma modes propagation properties at fre- 
quencies below the gap frequency and such properties are quite distinct below and above 
the c-axis plasma frequency. 

This work was done under a CNRS(France)-CNPq(Brasil) coUabotation program. 
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VII. APPENDIX 

In this appendix we provide further details on how the exact parametrized solutions of 
section |TV| leads the the pictures developed in sections fT^ and section |V| for below and 
above the plasma frequency, respectively. 

The three possible regimes, optical coupled and assymptotic follow from the exact 
parametrized solution. In particular we show that Eq.(p5D and Eq. (|26|) , for the coupled 



regime, and Eq. (p9D , for the asymptotic regime, follow from the Eq.(^) and Eq.(^). We 
derive both Eq. (^^ and Eq. (^) within the following approximations: (i) retardation effects 
are neglected (7^1) (VI — 7^ ~ 1); (H) fields evanesce slowly inside the film, t ^ 1. 
Therefore the Eq. ([T8|) dispersion relations become a/A X^y'^ ~ 1/2 and VA X^y'^ ~ 2/t^ for 
the symmetric and for the antisymmetric modes, respectively; and (in) the film thickness 
is much smaller than the London penetration length along the surfaces A ^ 1 so that 



even though t is small we can approximate Eq.(1T7|) hj At'^ ^ JX R/{1 — 7^Xi?/r). Di- 
rect elimination of the parameter t leads to both Eq.(p5D and Eq. (pB]) . In the asymptotic 
regime fields inside the superconductor fade away very quickly from the surfaces , thus cor- 
responding to t ^ 1. The surface decoupling is seen from the antisymmetric and symmetric 
dispersion relations which converge to the same asymptotic frequency ujas- According to 
Eq.(pi?D gA{t —>■ 00) = gs[t -^ cxd) — > t^ and Eq.(pPD follows in a straighforward way from 
this argument. 

To derive Eq. ( PPJ ) from the exact parametrized solution is in fact very simple, because it 
does not demand a detailed study of the 7(t') curve. To do so we reparametrize the exact 
solution of Eq . (|23|) , replacing t' by a new parameter ti, t' = tj — ti where ts = {2 N + l)7r, 
and t^ = 2 (A^ + 1) vr. The advantages are two-fold: both symmetric and antisymmetric 
modes are now defined in the same interval < ti < vr; and the approximation is under 
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control, corresponding to the limit ti/ti ^ 1. Notice that at ti = all curves satisfy 
the condition 7 = 1, thus justifying our claim that Eq.(pOD is a good approximation for 
the exact parametrized dispersion relation of Eq.(p3D and Eq.(p^ when the phase velocity 
is approximately given by the speed of light in the dielectric. Indeed to obtain Eq.(pOD 
from Eq.([23|) and Eq.(P^, notice that 7(t7 — ti) can be Taylor expanded in powers of ti 
because the functions hj(t'), introduced in section |^ are well behaved in this neighborhood: 
hj{t') = (ti — ti)^tan^ (ti/2). We introduce the approximation ti/ti < 1 into Eq.(p4D 
obtaining that X/ ^ {r/R)/{'y] — r/[yl(t/)^ + 1)]}. The key issue here is that all the 
parametric dependence of X is now on 7, because the last term was approximated by 
a constant, r/A{ti — ti)^ ~ r/Atj. Therefore we obtain the approximated parametrized 



dispersion relation. 



Uj{ti) ^ Up 



li{ii, 



^i)^-z*j 



(33) 



- (34) 



'^1? 



from where the family of curves oj{(3) of Eq.(|30D follow, by suitably removing the function 
7(tl). 

Finally we would like to get some more detailed information about the parametrized exact 
solution, in particular the two following issues: (^zj ^^{{3) is always an increasing function of 
/?; (ii) near ujp it suffices to consider the positive root of Eq.(|23|) because the positive and 
negative roots of Eq.(P3[) are just parts of the same curve and meet where the square root 
vanishes. The claims are of general validity but we derive them under some further working 
assumptions on Eq.(E3), that in leading order, acquires a very simple form. 



1 ± y'l - 45J tan2 (I) 



Ait,-t,r>.r»l, ,K«^ = ^_^^_|^^_^ (35) 



To have the condition r ^ 1 satisfied one must choose a dielectric medium of sufficiently 
high constant: e ^ Eg- To have at least one mode described by the above inequality, one 



must require that the minimal frequency satisfies the inequality: t^^^ ^ v ^M? which is 
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eJEs ^ 27rA||/a. To show that the curve never saturates, take t\ = 0, where one finds 
that 7=1 and 7 = for the positive and negative roots, respectively. Obviously the 
positive root is the low frequency part of the curve. Notice that the negative root is not 
physical because for 7 = P{X) is negative according to Eq. (|2^) . This just means that 
the negative square root curve must end at a finite non-zero 7, where the denominator of 
Eq.(^4|) vanishes and so X diverges. 7 is limited between zero and one, so it must also be 
finite at this point. Consequently u must also diverge, and so, we conclude from this that 
the curve never saturates in a; or /5. 
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FIGURES 
FIG. 1. A pictorial view of wave propagation in a superconducting filni surrounded by two 

equivalent non-conducting media. For oj < ujp the instantaneous electric field and the superficial 

charges are shown for symmetric (a) and antisymmetric (b) modes. For uj > u)p the optical ray 

associated to the plane wave travelling inside the film is shown here (c) for M=l,2,3 and 4. The 

symmetry of each state is also shown here, for both cases (a), (b) and (c), through the diagram Ez 

versus x. 

FIG. 2. The symmetric and antisymmetric uj < ujp dispersion relations are shown for three 
anisotropics. A very thin superconducting film, d = 10n?7i-thick, surrounded by SrTiO^, is con- 
sidered. 

FIG. 3. This Figure shows the asymptotic frequency versus e, for three anisotropics. For small 
e the asymptotic frequency is just cOp. The asymptotic frequency for YBCO is the most sensitive 
to changes in the susbstrate dielectric constant. 

FIG. 4. The dispersion relations above Wp is shown for a 100n?7i-thick film and anisotropy 
'^±/'^|| = 100. The M = 1 mode propagates with speed v and at the M = 66 mode is the upper 
limit for the validity of the present theory. 
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